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I.  INTRODUCTION 


The  SOLEX  solar  X-ray  experiment  aboard  the  U.  S.  Air  Force  Space 
Test  Program  P78-1  satellite  consists  of  two  collimated  Bragg  crystal 
X-ray  spectrometers  (Landecker,  McKenzie,  and  Rugge  1979;  McKenzie  et 
al.  1980b).  The  experiment  was  designed  to  study  solar  activity, 
primarily  through  observations  of  individual  active  regions.  Measure¬ 
ments  of  solar  flare  spectra  have  been  hampered  by  the  combination  of 
small  spectrometer  fields  of  view  (20  and  60  arc  sec  FWHM)  and  the 
lack  of  a  real  time  pointing  capability.  Nevertheless  we  have  ob¬ 
tained  spectra  over  the  7.8  -  23.0  A  range  during  the  rise  phase  of  a 
flare  on  1979  June  10  (McKenzie  al.  1980a, b)  and  spectra  of  the 
density-sensitive  0  VII  lines  for  flares  on  1980  April  8  and  May  9 
(Doschek  et  al.  1981),  The  present  paper  is  an  analysis  of  a  series 
of  13  X-ray  spectra  in  the  7.8  -  23.0  A  range  taken  primarily  during 
the  decay  phase  of  a  flare  on  1979  March  31. 

A  time  series  of  X-ray  spectra  provides  the  opportunity  to  ana¬ 
lyze  the  evolution  of  the  hot  coronal  plasma  during  a  flare.  In  con¬ 
trast  to  the  situation  for  the  1979  June  10  flare  (McKenzie  _et_  al. 
1980b),  for  which  we  obtained  spectra  only  during  a  relatively  long 
rise  phase,  for  the  1979  March  31  flare  we  have  spectra  covering  ap¬ 
proximately  2100  seconds  during  part  of  the  rise,  the  maximum,  and 
decay  phase  of  the  event.  The  X-ray  emission  peaked  during  the  second 
spectral  scan.  While  the  June  10  flare  data  showed  no  major  changes 
in  the  temperature  structure  of  the  plasma  during  the  rise  phase,  the 


March  31  data,  as  expected,  show  the  flare  cooling  during  the  decay 
phase. 


II.  THE  FLARE  DATA 


We  consider  X-ray  observations  of  McMath  Plage  Region  15918 
(active  region  1661)  on  1979  March  31  between  2321  and  2356  UT.  At 
the  start  of  our  observations  a  flare  was  in  progress.  According  to 
data  returned  by  the  MONEX  LEM  (low  energy  monitor),  a  full  disk  view¬ 
ing  proportional  counter  operating  in  the  range  1-22  keV  (0.6  - 
12.4  A)  also  aboard  P78-1  (Landecker,  McKenzie,  and  Rugge  1979),  the 
peak  X-ray  emission  occurred  between  2323.2  and  2324.5  UT  with  higher 
energy  channels  peaking  earlier.  The  SOLRAD  11  1-8  A  detector 
recorded  a  maximum  flux  of  about  3  x  10~^  erg-cm^-s-*  (X-ray  class 
M3),  according  to  plots  published  in  Solar  Geophysical  Data  (#417, 
Part  I).  An  Ha  flare  observed  in  McMath  Region  15918  at  coordinates 
S24,  E19  started  at  about  2318  UT,  peaked  at  2321,  and  ended  after 
2400  UT,  and  was  assigned  importance  IB  by  observers  at  one  observa¬ 
tory  and  2B  at  three  others  (Solar  Geophysical  Data  #416,  Part  I). 
The  P78-1  tape-recorded  observations  ended  at  2356.2  UT,  by  which  time 
the  X-ray  emission  had  decayed  by  at  least  65%  from  the  peak  in  all  of 
the  spectral  lines  of  highly-ionized  iron  to  be  considered  below. 

The  SOLEX  spectrometers  were  scanning  at  a  rate  of  31.25  30.2- 
arc-second  steps  per  second  between  Bragg  angles  of  17. 4®  and  61.7° 
during  the  March  31  flare.  SOLEX  B,  collimated  to  60”  x  60"  (FWHM), 
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used  a  rubidium  acid  phthalate  crystal  (RAP;  2d  ■  26*12  A)  to  cover 
the  7.8  -  23.0  A  range.  SOLEX  A,  collimated  to  20M  x  20"  (FWHM),  used 
ammonium  dihydrogen  phosphate  (ADP;  2d  *  10.64  A)  to  cover  the  range 
3.2  -  9.4  A.  The  time  required  to  record  a  complete  spectrum  was 
168.96  s,  and  each  line  was  sampled  only  very  briefly  during  this 
period.  Upon  reaching  one  end  of  a  spectral  jean,  the  spectrometer 
changed  direction  and  scanned  toward  the  other  programmed  end  point. 
Therefore,  the  time  interval  between  successive  samplings  of  a  given 
line  was  not  constant.  Lines  near  the  scan  ends,  such  as  the  0  VII 
lines  considered  below,  were  sampled  twice  in  rapid  succession  and 
then  were  unobserved  for  a  relatively  long  period  of  time.  In  addi¬ 
tion,  SOLEX  shares  the  pointing  section  of  the  OSO-like  P78-1  satel¬ 
lite  with  the  NRL  SOLWIND  white  light  coronagraph  (Sheeley  ^t_  al. 
1980).  At  intervals  of  about  ten  minutes  the  solar  experiments  are 
automatically  pointed  at  sun  center  so  that  a  SOLWIND  exposure  of 
about  one  minute  duration  can  be  made.  During  the  exposure  the  SOLEX 
scan  is  not  interrupted,  and  the  result  is  that  about  one-fourth  of 
the  spectra  have  gaps  in  them.  Because  of  the  small  field  of  view  of 
SOLEX  A  and  a  pointing  offset  of  28  arc  seconds  between  the  two  spec¬ 
trometers,  the  SOLEX  A  counting  rates  were  small;  hence  this  spectro¬ 
meter  provided  little  data  to  supplement  the  SOLEX  B  observations. 
Consequently  we  analyze  only  the  SOLEX  B  spectra  here. 

Figure  1  shows  selected  lines,  excited  at  different  temperatures, 
from  three  spectral  scans  during  the  decay  of  the  1979  March  31 
flare.  Data  are  corrected  for  detector  electronics  dead  time  but  not 


11 


COUNTS/C  032  S 


Figure  1:  Dead  time  corrected  profiles  of  selected  X-ray  lines 

observed  during  the  1979  March  31  flare*  Each  panel 
displays  measurements  from  one  spectral  scan  with  the 
time  at  which  each  line  was  observed  shown  above  the 
profile.  The  horizontal  bar  beneath  each  profile  is  at 
the  local  background  counting  rate  level.  The  lines  are 
identified  in  Table  1* 
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Instrument  response,  and  the  horizontal  line  beneath  each  profile  is 
at  the  background  level*  The  lines  and  the  approximate  temperatures 
at  which  they  are  most  efficiently  produced  are  shown  in  Table  1*  The 
plots  are  scaled  so  that  the  highest  32  ms  count  in  the  Fe  XVII  line 
is  nearly  full  scale.  The  more  rapid  decay  of  the  Fe  lines  from 
higher  stages  of  ionization  makes  apparent  the  cooling  of  the  flare 
plasma* 


III.  ANALYSIS 

The  spectral  lines  identified  in  connection  with  Figure  1  span  a 
wide  range  of  excitation  temperatures,  and  their  time-dependent  fluxes 
can  be  used  to  describe  the  cooling  of  the  flare.  In  Figure  2  the 
fluxes  of  a  selection  of  these  lines,  normalized  to  the  maximum  ob¬ 
served  flux  for  each  line,  are  plotted  as  a  function  of  time.  The 
sample  error  bars  reflect  counting  statistics  only.  Systematic  uncer¬ 
tainties,  especially  in  estimating  the  background,  could  increase  the 
error  bars  substantially.  Each  line  flux  was  computed  by  integrating 
dead  time  corrected  counts  under  the  line  profile  and  subtracting 
background.  Since  each  flux  is  compared  only  with  the  same  line's 
flux  at  a  different  time,  no  correction  for  instrument  response  was 
needed  in  constructing  Figure  2.  McKenzie  et  al.  (1980b)  give  a  de¬ 
tailed  description  of  this  data  analysis  procedure. 

Since  the  Fe  XIX  lines  plotted  in  the  figure  are  blended  with  the 
Ne  IX  ls^  lS0  “  ls2p  line,  a  correction  had  to  be  made.  This  was 
done  by  using  the  flux  of  the  Ne  IX  ls^  *Sq  -  Is 2s  "*Sj  line  at 
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TABLE  1 


Selected  X-ray  Linas 


Transition  Wavelength  Temperature  Remarks 


(X)  uo6k> 


o  vu 

l.2  1S0  -  1.2p  lPx 

21.60 

2 

o  vm 

X*  ZS1/2  -  2p  2P1/2>3/2 

18.  97 

3 

Fe  XVU 

2p6  1SQ  -  2p53d  lPt 

15.01 

4 

Fe  XVIH 

P3/2  "  ^  D)3d  P3/2 

14.21 

6 

blend 

Fe  XIX 

2p4  3P2  -  2p3(2D)3d  3Dj,  3P2 

13.52 

7 

blend,  also  with 

Ne  DC  l*2  1S0  - 
l«2p  3P,  (•«•  tea) 

Fe  XX 

2p3  4S3/2-2p2(3P)3d4P5/2i3/2  l/2 

12.83 

9 

blend 

Fe  XXI 

2p2  3P0  -  2p3d  3Dj 

12.28 

10 

blend  with  Fe  XVU 

2p6  lSQ  *  2p54d 

3_.  0 

Dl 

Fe  XXIV 

2i  2S1/2  -  3p  2P3/2>  1/2 

10.  64, 

10.66  16 

partial  blend, 
treated  as  one  line 
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FLUX /(MAXIMUM  FL 


Figure  2:  The  line  flux  of  selected  lines  of  highly  ionized  iron, 

normalized  so  that  the  maximum  observed  flux  is  1.00. 
The  measurements  are  connected  to  facilitate  reading  the 
plots.  The  lines  are  identified  in  Table  1. 


15 


13.701  A.  Gabriel  and  Jordan  (1969,  1972)  showed  that  the  flux  ratio 
of  the  to  the  line  is  given  by 


1.35  N  C(2  3S  -  2  3P) 


R  -  0.33  (l  +  - \ - = - 

A(2  JS  -  1  XS) 


where  Ng  is  the  electron  density.  By  taking  a  recently  calculated 
value  for  the  collisional  excitation  coefficient  (Pradhan,  Norcross, 
and  Hummer  1981)  and  A(2  -  1  *S)  tabulated  by  Gabriel  and  Jordan 

(1972),  we  find,  for  a  temperature  of  4  x  10^  K, 


1  +  1.30  x  10_12  N 


We  do  not  believe  that  R  for  neon  varied  significantly  from  its  low 
density  limit  of  3.0  during  this  flare.  We  show  below  that  the  maxi* 
mum  observed  density  for  the  0  VII  emitting  plasma  at  2  x  10^  K  was 
4.7  x  10^  cm”**.  Equation  (2)  shows  that  a  density  of  7.7  x  10^  cm*'* 
at  4  x  10^  K  would  be  required  for  the  derived  Ne  IX  ^P^  flux  to  be 
erroneous  by  10%.  Although  such  errors  cannot  be  ruled  out  complete¬ 
ly,  the  corrections  to  the  Fe  XIX  fluxes  were  small,  ranging  from  less 
than  10%  before  2342  UT  to  about  25%  for  the  last  plotted  measure¬ 
ments.  Thus  the  Fe  XIX  decay  curve  would  not  be  significantly 
affected  by  errors  in  the  Ne  IX  JP^  flux. 

From  Figure  2  one  sees  immediately  that  the  fluxes  of  all  of  the 
plotted  lines  began  to  decay  within  a  few  minutes  of  the  X-ray  flux 
maximum.  Lines  from  0  VII,  0  VIII,  and  Fe  XXI,  not  plotted  to  avoid 
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severe  crowding  in  the  figure,  show  similar  decay  curves.  0  VII  lines 
are  emitted  at  2  x  10^  X,  and  Fe  XXIV  lines  at  about  1.6  x  10^K.  Work 
by  Feldman  et  al.  (1980;  see  also  Doschek  et  al.  1980)  has  shown  that, 
at  least  after  2326  UT,  very  little  of  the  plasma  in  this  flare  was  at 
temperatures  above  1.8  x  10^  K.  Thus  the  fluxes  of  lines  emitted  at 
all  temperatures  present  in  the  coronal  flare  began  to  decrease  a 
short  time  after  the  time  of  maximum  X-ray  flux. 

The  simultaneous  decay  of  all  the  observed  emission  line  fluxes 

indicates  that  the  flare  emission  measure,  the  volume  integral  over 

2 

the  flare  of  the  square  of  the  electron  density  (/  dV),  also  began 
to  decrease  a  short  time  after  flare  maximum.  This  is  contrary  to 
results  based  upon  observations  made  with  broad  band  detectors  (Horan 
1971)  which  have  the  emission  measure  continuing  to  increase  long 
after  the  peak  X-ray  flux  is  observed.  If  we  take  the  time  of  the  X- 
ray  peak  to  be  2323  UT,  at  which  time  the  6  keV  X-ray  flux  measured  by 
the  MONEX  LEM  reached  a  maximum,  observations  require  that  the  emis¬ 
sion  measure  began  to  decrease  no  more  than  500  s  after  the  X-ray 
peak.  Doschek  et  al.  (1980)  observed  that  the  emission  measure  for  a 
number  of  large  flares  peaked  only  a  short  time  after  the  X-ray  flux 
maximum.  Their  observations  were  of  lines  of  Ca  XIX  and  Fe  XXV,  spe¬ 
cies  present  at  temperatures  above  10^  K.  The  present  observations 
extend  the  result  over  the  entire  regime  of  temperatures  in  the  coro¬ 
nal  flare  plasma.  The  apparent  error  in  the  broad  band  results  prob¬ 
ably  lies  in  the  method  of  analysis.  A  temperature  is  deduced  from 
the  broad  band  spectrum  and  then  emission  measure  is  determined  by 
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using  this  temperature  and  the  absolute  flux  in  a  selected  energy 
band.  The  deduction  of  temperature  is  by  no  means  easy,  and,  since 
fluxes  are  often  measured  at  energies  substantially  higher  than  kT,  a 
small  error  in  temperature  can  lead  to  a  large  one  in  emission  meas¬ 
ure.  Effects  such  as  pulse  pile-up  (Datlowe  1975a)  may  further  com¬ 
plicate  the  process.  The  present  result  is  unequivocal:  all  observed 
X-ray  line  fluxes  decreased  simultaneously  so  the  coronal  emission 
measure  must  have  decreased. 


We  define  the  differential  emission  measure,  e(T),  by  considering 
the  volume  of  flare  material  emitting  at  temperature  T.  The  total 

0  AV 

emission  measure  can  be  written  as  /  ^  dT.  We  define  e(T)  as 

2  dV 

Ne  dx*  use  the  measured  line  fluxes  to  compute  e(T)  for  the 

coronal  flare  plasma.  To  establish  the  correctness  of  the  results  we 
will  describe  in  detail  the  atomic  data  used  and  some  inflight  instru¬ 
ment  calibration  checks. 


The  volume  emission  rate  for  the  spectral  line  j  ♦  i  at  tempera¬ 
ture  T,  arising  from  collisional  excitation  from  the  ground  state,  is 


P 


ji 


8.63xl(T6  N  N  8  .  B..  hv  exp(-E  VkT) 

_ e  6  8J  J1 _ _ _ SJL _ _ 

u  /T 
g 


-3 

erg-cm 


-1 

-s 


(3) 


where  Ng  is  the  ground  state  population  (cm  )  of  the  emitting 
ion,  ftgj  is  the  collision  strength  from  the  ground  state  to  state  j 
(averaged  over  a  Maxwellian  electron  distribution),  is  the  branch¬ 
ing  ratio,  hv  the  X-ray  energy,  Egj  the  energy  of  state  j,  and  the 
statistical  weight,  2Jg+l.  If  the  differential  emission  measure  and 
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the  collision  strengths  are  known,  the  line  fluxes  may  be  computed 
from  equation  (3).  Conversely,  if  enough  line  fluxes  and  collision 
strengths  are  known,  the  differential  emission  measure  can  be  esti¬ 
mated. 

The  difficulty  in  using  equation  (3)  lies  in  the  scarcity  of 
collision  strength  data.  Collision  strengths  are  difficult  to  calcu¬ 
late,  and,  in  the  X-ray  region,  few  are  known.  Workers  at  Los  Alamos 
Scientific  Laboratory  have  compiled  calculated  collision  strength  data 
from  a  large  number  of  sources  (Magee  et_  al.  1977;  Merts  et_  al. 
1980).  The  results  of  calculations  of  collision  strengths  between  LS 
terms  are  given  in  the  form  of  a  convenient  fitting  formula.  Colli¬ 
sion  strengths  for  individual  J  levels  are  not  given.  We  have  used 
the  compiled  results  for  the  transitions  of  0  VII,  0  VIII,  Fe  XVII,  Fe 
XVIII,  Fe  XXI,  and  Fe  XXIV  shown  in  Figure  1.  Collision  strengths  for 
the  other  lines  in  Table  1  are  unavailable.  We  discuss  the  use  of 
these  data  in  the  following  paragraphs. 

The  use  of  the  collision  strength  data  is  straightforward  in  the 
cases  of  0  VII,  0  VIII,  Fe  XVII,  and  Fe  XXIV.  For  0  VII  and  Fe  XVII 
we  treat  singlet-singlet  transitions  so  the  Los  Alamos  compiled 
results  can  be  used  directly.  In  addition  to  collisional  excitation 
from  the  ground  state,  cascade  from  the  2s2p^3d  level  populates 
the  Fe  XVII  2s^2p^3d  level  (Loulergue  and  Nussbaumer  1973).  The 
collision  strength  for  the  *D2  level  was  taken  from  Merts  et_  al. 
(1980).  For  0  VIII  the  two  transitions  are  completely  unresolved,  so 
the  entire  collision  strength  is  used.  For  Fe  XXIV  the  two  lines  are 
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partially  blended ,  and  the  numerical  analysis  is  more  easily  accom¬ 
plished  if  they  are  treated  together.  The  collision  strength  is  in 
excellent  agreement  with  recent  calculations  by  Hayes  (1979). 

For  Fe  XXI  the  situation  is  complicated  by  the  fact  that  the 
3  3 

upper  term  is  D  and  the  ground  term  is  (the  ground  state  is 

2p  Pq)«  There  are  six  dipole  allowed  transitions  between  the  two 
multiplets.  Mason  et  al.  (1979)  have  calculated  collision  strengths 
for  Fe  XXI  individual  J  levels  at  energies  of  80,  100,  and  150 

Rydbergs.  Their  summed  LS  term  to  term  collision  strengths  are  in 
excellent  agreement  with  those  of  Merts  et  al.  (1980).  We  have  used 
the  Merts  _e£  al.  results  because  of  their  convenient  form  and  have 

apportioned  the  LS  collision  strengths  among  the  J  transitions  in 
proportion  to  the  Mason  £t^  al.  80  Rydberg  (1.09  keV)  results.  80 
Rydbergs  is  near  the  threshold  energy  for  the  transitions.  The 
branching  ratio  for  the  2p2  **Pq  -  2p3d  transition,  0.845,  was 

determined  from  results  in  Mason  e^  al.  Data  showing  this  line  and 
the  competing  2p2  ^Pj  -  2p3d  transition  at  12.40  A  (McKenzie  et 

al. ,  1980b)  verify  that  this  branching  ratio  is  approximately  correct. 

The  Fe  XXI  line  flux  had  to  be  corrected  for  the  contribution  of 
the  partially  blended  Fe  XVII  2p^  *Sq  -  2p^4d  line  (12.264  A; 

Hutcheon,  Pye,  and  Evans  1976).  Figure  1  clearly  shows  that  this 
blend  was  quite  Important  late  in  the  flare.  The  correction  was  made 

by  determining  the  4d  Dj  line  flux  late  in  the  flare,  when  it  domi¬ 

nated  the  Fe  XXI  line,  and  comparing  it  to  the  Fe  XVII  2p^  *Sq  - 
2p^3d  flux  at  that  time.  The  ratio  of  fluxes  of  the  two  Fe  XVII 
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lines  (0*086)  was  then  used  to  correct  the  earlier  data*  Loulergue 
and  Nussbaumer  (1973,  1975)  have  shown  that  this  line  ratio  is  only 
weakly  dependent  on  T  provided  that  T  >_  5  x  10^  K. 

The  most  difficult  problem  was  with  the  Fe  XVIII  blend,  2p5  ^3/2 
-  2p^(*D)3d  ^5/2*  ^3/2  (Feldman  et  al*  1973).  The  Los  Alamos  compi¬ 
lation  includes  LS  collision  strengths  for  both  transitions  but  we 
know  of  no  such  calculations  for  individual  J  levels.  We  therefore 
repotted  to  the  approximation. 


8tt 
✓  3 


fij8  » 


(4) 


where  f^j  is  the  absorption  oscillator  strength  and  g  the  Gaunt  fac¬ 
tor.  Thus,  for  two  transitions  having  the  same  and  g, 


ai  (2Jj  +  1)  sl 

“  (2J2  T  l)  f2  * 


where  the  J's  refer  to  the  lower  state  of  each  transition.  The  Fe 
XVIII  collision  strengths  were  approximated  by  using  the  LS  collision 
strengths  from  the  Los  Alamos  compilation  and  equation  (5)  with  oscil¬ 
lator  strengths  computed  by  Cowan  (1973).  The  branching  ratios  were 
derived  from  transition  probabilities  computed  by  Cowan. 

The  only  provision  for  inflight  calibration  of  the  SOLEX  spectro¬ 
meters  is  a  radioactive  source  that  uses  alpha  particles  to  generate 
Mg  Ka  X-rays  (9.9  A)  by  fluorescence.  From  time  to  time  the  spectro¬ 
meters  are  placed  in  a  special  configuration  so  that  the  detectors  can 
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be  irradiated.  These  calibration  runs  have  shown  that  the  SOLEX  B 
sensitivity  at  9.9  A  has  not  significantly  changed  since  before 
launch.  In  addition,  the  coronal  X-ray  spectrum  includes  a  number  of 
line  pairs  having  flux  ratios  that  are  only  weakly  dependent  on  tem¬ 
perature.  By  studying  these  line  pairs  in  active  region  (nonflaring) 
spectra,  we  have  been  able  to  verify  the  SOLEX  B  calibration. 

We  have  analyzed  line  pairs  from  five  SOLEX  B  active  region  spec¬ 
tra  recorded  in  March  and  April  1979.  The  flux  ratio  of  the  two  0  VII 
lines,  ls^  *Sq  -  ls3p  and  ls^  *Sq  -  ls2p  at  18.63  A  and 
21.60  A,  respectively,  depends  weakly  on  temperature.  We  find  an 
average  photon  flux  ratio  of  0.120  ±  0.022.  Rocket  flight 
measurements  of  an  active  region  by  Parkinson  (1975)  and  the  quiet 
corona  by  McKenzie  ^t^  al.  (1978)  found  photon  ratios  of  0.129  and 
0.103,  respectively.  Acton,  Catura,  and  Joki  (1975)  pointed  out  that 
the  flux  ratio  of  the  0  VIII  Is  ^$1/2  ~  2p  ^pi/2  3/2  l*ne 
(18.97  A)  to  the  He  IX  Is2  -  ls2p  1P1  (13.45  A)  line  is  only 
weakly  temperature  dependent.  From  25  rocket  observations  they 
obtained  a  photon  flux  ratio  broadly  distributed  with  an  average 
between  7  and  8.  The  five  SOLEX  spectra  yield  an  average  ratio  of  6.1 
±  0.9.  Parkinson  (1975)  found  6.8.  The  ratio  of  the  Fe  XVII  lines, 
2p^  *Sq  -  2p^3s  *Pi  at  16.78  A,  and  2p^  *Sq  -  2p^3d  *P^  at  15.01  A  is 
very  weakly  temperature  dependent,  according  to  Loulergue  and 
Nussbaumer  (1975).  They  cite  data  from  eight  different  spectra.  The 
simple  average  of  the  energy  flux  ratio  for  these  measurements  is  0.68 
±  0.14.  The  five  SOLEX  active  region  spectra  yield  an  average  energy 
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flux  ratio  of  0.69  ±  0.02.  The  Fe  XVII  ratio  2p6  -  2p54d  3DL 

(12.26  A)  to  2p^  *Sg  -  2p33d  (15.01  A)  is  more  temperature 

dependent  than  the  other  ratios  cited  above,  but  the  dependence  is 
weak  at  temperatures  for  which  Fe  XVII  lines  are  strongly  emitted. 
The  average  photon  flux  ratio  for  the  SOLEX  B  measurements  was  0.065  ± 
0.005,  in  good  agreement  with  0.065  of  Parkinson  (1975),  but  well 
above  the  values  of  Hutcheon,  Pye,  and  Evans  (1976).  Walker,  Rugge, 

and  Weiss  (1974)  find  a  ratio  of  0.061,  but  they  caution  that  the  4d 

line  may  be  blended  in  their  spectra.  Taken  together,  the 
measurements  cited  above  verify  the  SOLEX  B  calibration  over  a  wave¬ 
length  range  including  all  of  the  lines  studied  in  this  paper  except 
for  the  Fe  XXIV  lines  at  around  10.65  A.  Since  both  the  detector 
efficiency  and  the  RAP  reflectivity  (Burek  1976)  are  uncomplicated 
down  to  the  Mg  K  absorption  edge  at  9.5  A  (the  detector  has  a  MgF2 
photocathode)  we  are  confident  in  the  calibration  above  this  wave¬ 
length. 

The  differential  emission  measure  was  computed  with  the  iterative 
procedure  of  Sylwester,  Schrijver,  and  Mewe  (1980).  In  addition  to 
the  collision  strengths  discussed  above,  we  used  ionization  equilibri¬ 
um  calculations  by  Jacobs  al.  (1977,  1978)  and  the  coronal  abun- 
dances  (A(Fe)  -  5.25  X  10-5,  A(0)  -  4.27  x  10“4,  A(H)  -  1)  of  Withbroe 
(1976).  e(T)  was  constrained  to  be  zero  at  and  above  1.8  x  10^  K,  in 
accord  with  the  results  of  Feldman  jat_  al.  (1980)  for  this  flare. 
Figure  3  shows  the  results.  Each  curve  must  be  divided  by  fc,  the 
collimator  transmission  averaged  over  the  flaring  region.  fc  ■  1  for 
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a  point  source  on  axis,  and  iQ  ■  0*25  for  a  flare  uniformly  filling 

the  collimator  field  of  view.  Smaller  values  occur  if  the  flare  is 

well  off-axis.  Examination  of  a  raster  map  of  the  active  region  made 

82  minutes  after  the  flare  peak  indicates  that  the  strongest  emission 

was  about  10  arc  sec  off  axis  at  that  time.  If  we  assume  that  this 

point  corresponds  to  the  center  of  the  flare  and  correct  for  solar 

rotation,  then  the  collimator  was  offset  by  about  20  arc  sec  for  the 

flare  observations.  Based  upon  this  offset  and  the  assumption  that 

the  characteristic  size  of  the  flare  is  unlikely  to  have  exceeded  40 

arc  sec  (Landecker  and  McKenzie  1980),  we  estimate  that  f  *  was  in  the 

c 

range  1.50  -  1.75. 

At  temperatures  below  4  x  10^  K,  e(T)  is  uncertain  because  of 
uncertainties  in  the  0/Fe  abundance  ratio.  Because  the  0  VIII  line 
has  significant  emissivity  at  temperatures  as  high  as  10^  K,  the  data 
summarized  in  Figure  3  can  be  used  to  estimate  the  0/Fe  abundance 
ratio.  We  were  able  to  obtain  satisfactory  fits  ,  fits  that 

reproduced  all  line  fluxes  to  within  15%)  for  A(0)/A(Fe)  in  the  range 
5  -  10.  Withbroe  (1976)  lists  values  of  this  ratio  for  the  corona 
ranging  from  4.6  to  27  and,  in  his  final  tabulation  of  coronal  abun¬ 
dances,  A  (0)/A(Fe)  *  8.1.  Figure  4  illustrates  the  uncertainty 

in  e(T)  for  one  of  the  spectra.  We  have  plotted  the  range  of  e(T) 
values,  basing  the  estimated  error  on  the  "continuous”  error  case  of 
Sylwester,  Schrijver  and  Mewe  (1980).  Two  plots  are  shown,  for 
A(0)/A(Fe)  of  5  and  10.  The  uncertainty  introduced  by  lack  of 
knowledge  of  the  abundances  is  very  large  below  4  x  10^  K  and  does  not 
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An  illustration  of  the  effect  of  element  abundances  on 
the  derived  differential  emission  measure  function.  The 
range  of  acceptable  values  of  e(T)  at  2334.5  UT  is 
plotted  for  two  assumptions  regarding  the  O/Fe  abundance 
ratio.  As  a  result  of  the  uncertainty  in  this 
ratio,  e(T)  is  almost  unknown  for  T  below  4  x  iO^K. 


vanish  until  about  8  x  10^  K.  We  also  note  that  errors 
in  e(T)  confined  to  a  narrow  temperature  range  might  be  considerably 
larger  than  the  uncertainties  displayed  in  Figure  4  (Craig  and  Brown 
1976;  Sylwester,  Schrijver,  and  Mewe  1980). 

The  outstanding  features  of  the  curves  in  Figure  3  are  the 
large  e(T)  at  low  and  high  temperatures  and  the  small  e(T)  at  around 
107  K  that  develop  several  minutes  after  the  flare  peak.  The  low 
emission  measure  at  1G7  K  is  manifested  by  low  Fe  XXI  emission.  We 
have  verified  that  this  low  emission  does  not  result  from  overcorrect¬ 
ing  for  the  Fe  XVII  2p^  *Sq  -  2p^4d  line.  As  noted  above,  we 

measure-  u,  -*D|)/(2p^3d  ratio  of  0.086  late  in  this  flare 

and  a  ratio  of  0.063  in  a  number  of  active  region  spectra.  It  is  not 
surprising  that  the  ratio  is  larger  in  a  flare.  We  have  found,  how¬ 
ever,  that  the  dip  in  e(T)  als*.  appears  if  a  ratio  of  0.065  is  used  to 
correct  the  flare  data.  With  regard  to  the  low  temperature  e(T),  re¬ 
call  the  large  uncertainty  shown  in  Figure  4.  We  cannot  say  wheth¬ 
er  e(T)  peaks  at  around  5  x  10^  K  or  is  high  at  all  temperatures  below 
this  level. 

Figure  3  suggests  that  e(T)  consists  of  separate  components  at 
high  and  low  temperatures.  The  material  at  107  K  early  in  the  flare 
may  be  part  of  the  low  temperature  component  which  cools  relatively 
quickly  to  reveal  a  high  temperature  part  with  emission  measure  in¬ 
creasing  sharply  with  increasing  temperature.  In  studies  of  the  ener¬ 
gy  balance  of  active  regions  and  the  transition  zone,  Jordan  (1976, 
1980)  cites  evidence  for  e(T)  of  the  form  aTb  for  T  up  to  some  cut-off 
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temperature.  The  hot  component  appears  to  be  consistent  with  this 
form.  This  would  mean  that  the  hot  component  has  very  little  material 
below  about  10^  K.  Although  the  cool  material  shows  little  change  in 
emission  measure  during  our  observations,  it  should  not  be  thought  of 
as  mere  background.  The  0  VII,  0  VIII,  and  Fe  XVII  fluxes  near  the 
peak  of  the  flare  are  approximately  10  times  those  in  nonflaring  ac¬ 
tive  regions. 

By  integrating  the  functions  in  Figure  3  for  T  above  2  x  10^  K  we 
can  determine  the  emission  measure  for  the  entire  coronal  flare.  Most 
previous  work  has  derived  only  an  apparent  emission  measure  based  on 
an  assumed  single  relatively  high  temperature  for  the  flare  (Datlowe 
1975b;  Feldman  et  al.  1980).  The  present  method,  however,  does  have 
the  disadvantage  that  the  result  must  be  divided  by  an  uncertain 
factor  fc*  The  integrated  emission  measures  are  given  in  Table  2.  We 
also  show  separately  the  integral  over  the  hot  material  below  10^  K. 

The  ratio,  R,  of  the  Is2  *Sq  -  ls2s  ^  to  the  Is2  ^Sq  -  ls2p  ^ 
lines  of  the  helium-like  ions  can  be  used  as  an  electron  density  diag¬ 
nostic  (Gabriel  and  Jordan  1969,  1972).  We  have  measurements  of  R  for 
0VII  and  Mg  XI.  From  a  previous  discussion  herein  one  can  see  that 
our  Ne  IX  measurements  are  unsuitable  for  this  purpose.  For  0  VII  we 
have  (similar  to  Equation  (2)  above), 

Ne  -  3.5  x  1010  -1)  (6) 
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TABLE  2 


Flap#  Parameters 


Time  <UT) 

Unit 

2326. 7 

2328.8 

2334.5 

2340. 1 

1.  8x1 07K 
lc  ^«<T)dT 

2xl06K 

(1048  cm'3)* 

+3.2 

15.2 

-1.5 

+3.9 

12.9 

-2. 1 

+2.7 

8.0 

-0.9 

+2.4 

7.0 

-1.0 

lxl  07K 
fc  ^«(T)dT 

(1048  cm'3)* 

+  3.  1 

8.7 

-1.3 

+3.3 

8.4 

-1.4 

+2.5 

6.0 

-0.7 

+2.2 

5.4 

-0.8 

2x1 06K 

Ne(2xl06K) 

(1010  cm'3) 

... 

4.  8+  0. 5 

2.8+  0.5 

1.0+  0.4 

Pe(2xl06K) 

(dyne-cm*^) 

— 

13.3+  1.4 

7.7+  1.4 

2.8  +  1.  1 

fc  PR 

(1026  erg-*'1)* 

+3.8 

8.6 

-1.3 

+3.6 

7.5 

-1.5 

+2.6 
4.8 
-0.  8 

+2.3 

4.2 

-0.8 

*  Uncertainty  based  on  C(T)  functions  for  extreme  values  of  A(Fe)/A(0)  (see 
Figure  4). 
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This  differs  from  the  expression  we  have  previously  used  (McKenzie  et 
al#  1980a)  because  here  we  use  new  collision  strength  calculations  by 
Pradhan,  Norcross,  and  Hummer  (1981).  Our  earlier  derived  densities 
should  be  reduced  by  a  factor  of  1.9  to  correspond  to  the  use  of  these 
collision  strengths.  Electron  densities,  Ne,  and  pressures,  Pe,  (P0  ■ 
NekT)  for  an  assumed  temperature  of  2  x  10°  K  are  given  in  Table  2. 
The  maximum  observed  density  of  4.8  x  10^  cm”^  is  lower  than  the  9  x 
10  cm  obtained  when  our  1979  June  10  observations  (McKenzie  et  al. 
1980a)  are  corrected  by  the  above-mentioned  factor  of  1.9,  and  is  much 
lower  than  the  10^  cm"^  found  for  flares  on  1980  April  8  and  May  9 
(Doschek  et  al.  1981).  In  the  March  31  flare  we  had  no  0  VII  obser¬ 
vations  until  after  the  peak,  and  densities  may  have  been  considerably 
higher  earlier.  The  observed  rapid  decrease  in  density  shows  that  the 
cooler  coronal  material  was  active  in  the  flare  process  despite  the 
slow  variation  in  e(T). 

For  Mg  XI  (T  -  6  x  106  K), 

Ne  -  1.6  x  1013  -  1)  (7) 

(Gabriel  and  Jordan  1972).  For  11  scans  the  mean  value  of  R  was  2.31 
with  a  standard  deviation  of  the  mean  of  0.08.  The  lowest  measured  R 
was  1.97  ±  0.28,  and  this  was  very  late  In  the  flare  (2354.8  UT).  We 
do  not  find  any  positive  density  determination  but  set  an  upper  limit 
of  9  x  10  cm  from  the  above  measurement  and  equation  (7). 
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Observations  from  Sky lab  (Vorpahl  et~  al.  1975;  Kahler,  Krieger, 

and  Vaiana  1975)  have  shown  that  the  hot  flare  plasma  is  located  in 

arches  or  loops.  A  pressure  equilibrium  can  establish  itself  in  the 

time  required  for  sound  to  travel  the  length  of  a  loop  (Moore  et  al. 

1980).  Then  a  short  time  after  the  flare  begins  to  cool  the  dif'feren- 

2  dV 

tial  emission  measure  can  be  expressed  as  the  product  (T)  -pjr;  that 
is,  the  temperature  can  be  expressed  as  a  function  of  position  along 
the  loop.  It  is  usually  assumed  (Moore  et  al.  1980)  that  the  tempera¬ 
ture  is  highest  at  the  top  of  the  loop  and  the  flare  cools  by  radia¬ 
tion  and  conduction  down  the  loop.  The  total  thermal  energy  in  the 


loop  can  be  written  as  ^  ? 

N  VdV 

E  -  j  k  /  (Ng  +  Nj)  TdV  -2.8  k  /  NgTdV  -  2.8lt  /  - 


2  A.'It'11 

2.8k2  /  — |-g— 


where  is  the  ion  density.  Equation  (8)  assumes  the  helium  abun¬ 
dance  is  20%  of  the  hydrogen  abundance  and  the  ion  temperature  equals 
the  electron  temperature. 

If  P0(T)  is  known  the  thermal  energy  can  be  computed.  Unfortu¬ 
nately  we  know  only  the  average  pressure  for  temperature  around  2  x 
106  K;  i.e.»  near  the  bottom  of  the  loops.  It  would  be  convenient  to 
assume  Pe  is  a  constant.  If  the  loops  are  in  hydrostatic  equilibrium, 
a  pressure  gradient  is  required  to  support  the  material  near  the 
tops.  If  the  weight  per  unit  area  of  the  loop  material  is  small  com¬ 
pared  to  the  pressure  near  the  feet,  constant  Pe  may  be  a  good  approx¬ 
imation.  The  weight  per  unit  area,  w,  is  given  by 
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w 


1.3  / 


"H  gs 


N  «  “  1.3  »„g > 


/  JL 

J  P  A 


P  A 
e 


»e2  S  « 


(9) 


where  is  the  proton  mass,  gs  is  the  gravitational  acceleration  at 
the  surface  of  the  sun,  A  is  the  cross-sectional  area?  of  the  loop(s), 
and  the  integral  is  over  T  greater  than  2  x  iO^K.  The  multiplying 
factor  1«3  arises  from  the  choice  of  helium  abundance  of  20%  of  the 
hydrogen  abundance.  Checking  the  above-stated  criterion  for  consist¬ 
ency  with  the  assumption  that  is  constant  we  require  w  «  Pe,  or 

A  >  1*3  /TNe2  *  (10) 

e 

provided  that  the  tubes  are  not  tapered •  For  the  March  31  flare  at 
2328.8  UT  this  becomes  A  »  5  x  1018  cm2  or  A  »  3.3  deg2.  The  Big 
Bear  Observatory  measured  the  Ha  flare  area  at  4.8  deg2  (Solar 
Geophysical  Data  #416,  Part  I).  If  we  take  the  Ha  area  as  an  upper 
limit  on  the  cross-sectional  area  of  the  X-ray  loops  we  see  that  the 
criterion  in  equation  (10)  is  not  satisfied;  the  weight  per  unit  area 
is  not  negligible  compared  to  the  pressure  near  the  base  of  the 
loops.  Since  T  enters  the  integral  of  equation  (8)  in  the  second 
power,  our  ignorance  of  the  pressure  at  high  temperature  is 
particularly  unfortunate.  Very  little  information  is  available  on  the 
high  temperature  flare  energetics.  Because  currently  known  density- 
sensitive  lines  excited  at  high  temperatures  only  enable  very  high 
densities  to  be  determined,  it  is  doubtful  that  this  situation  can  be 
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improved  in  the  near  future*  X-ray  or  XUV  images  with  very  high 
spatial  resolution  might  allow  densities  to  be  deduced,  but  such 


images  are  currently  unavailable* 

The  rate  of  radiative  energy  loss,  PR,  can  be  determined  easily: 

PR  -  /  A(T)Ng2  dT  erg/s  .  (11) 


A(T)f  the  radiative  energy  loss  rate  is  defined  by  the  equation. 
Table  2  gives  PR,  which  was  computed  using  A(T)  from  Raymond,  Cox,  and 
Smith  (1976),  corrected  for  the  higher  iron  abundance  we  use.  A  very 
rough  estimate  of  the  energy  content  can  be  made  by  performing  the 
integral  in  Equation  (8)  assuming  Pe  *  constant.  At  2328.8  UT  this 
gives  5  x  10  y  ergs  and  a  radiative  decay  time  of  7  x  10  s.  Only  if 
Pe  is  higher  at  high  T  (near  the  tops  of  the  loops)  than  at  low  T, 
will  the  radiative  decay  time  be  shorter.  In  this  case  hydrostatic 
equilibrium  would  not  obtain.  Moore  et  al.  (1980)  give  a  simple  ex¬ 
pression  for  the  conductive  decay  time: 


T 

C 


4  x  10 


-10 


s. 


(12) 


If  L  ■  2  x  109  cm,  a  typical  value  reported  by  Vorpahl  et  al.  (1975), 
Ne  *  5  x  1010  cm”^,  and  T  *  1.8  x  10^  K,  -  58  s.  Moore  et  al. 
estimate  that  might  be  an  order  of  magnitude  larger  if  the  flux 
tubes  are  tapered.  Although  the  estimates  are  quite  rough,  it  appears 
that  conductive  cooling  played  a  more  important  role  than  radiative 
cooling  in  this  flare. 


IV.  SUMMARY 


We  have  obtained  13  X-ray  spectra  (7.8  -  23.0  A)  during  a  solar 
flare  on  1979  March  31.  The  temperature  evolution  of  the  hot  flare 
plasma  can  be  seen  by  examination  of  this  time  series  of  spectra.  The 
differential  emission  measure  function  has  been  derived  at  four  times 
during  the  flare  decay.  The  conclusions  of  the  study  are  summarized 
in  this  section. 

The  fluxes  of  spectral  lines  emitted  at  temperatures  from  2  x  10b 
K  to  18  x  10bK,  the  maximum  flare  temperature,  decreased  within  500  s 
of  the  time  of  peak  X-ray  emission.  This  demonstrates  that  the  total 
emission  measure  also  began  to  decline  near  the  time  of  flare  maximum, 
contrary  to  results  derived  from  earlier  broad  band  observations.  The 
differential  emission  measure  had  two  components,  above  and  below 
about  107  K.  The  low  temperature  part  had  the  larger  emission  measure 
throughout  the  flare  decay  and  dropped  off  rapidly  at  temperatures 
above  5  x  10b  K.  The  high  temperature  component  rose  steeply  with 
temperature  in  the  range  10-18  x  10b  K  and  was  consistent  with  the 
form  e(T)  -  aTb. 

Knowledge  of  the  differential  emission  measure  function  enabled 
us  to  determine  the  total  emission  measure  and  the  radiative  energy 
loss  rate  from  the  flare.  The  electron  pressure  (or  density)  could  be 
measured  only  at  a  temperature  of  2  x  10b  K  where  a  ratio  of  0  VII 
lines  is  useful.  There  is  currently  no  satisfactory  way  to  determine 
densities  at  temperatures  above  a  few  million  K.  Furthermore,  in  this 
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flare  at  least,  there  is  no  justification  for  the  use  of  the  simplify¬ 
ing  assumption,  Pe  ■  constant.  A  consequence  is  that  very  little 
information  is  available  on  the  flare  energetics.  Despite  large  un¬ 
certainties  arising  from  our  ignorance  of  Pe  and  of  the  flare  geomet¬ 
ry,  we  conclude  that  conductive  cooling  was  dominant  over  radiative 
cooling  in  this  flare. 
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laboratory  operations 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmoa- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  g.aphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetoepheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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